Entanglement Transfer from Bosonic Systems to Qubits 
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We study the entanglement of a pair of qubits resulting from their interaction with a bosonic 
system. Here we restrict our discussion to the case where the set of operators acting on different 
qubits commute. A special class of interactions inducing entanglement in an initially separable two 
qubit system is discussed. Our results apply to the case where the initial state of the bosonic system 
is represented by a statistical mixture of states with fixed particle number. 
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INTRODUCTION 



The quantum correlations between subsystems present in entangled states are indispensable for many quantum 
communication protocols However, these correlations cannot be created by local operations and classical commu- 
nication(LOCC). Therefore, in order to entangle two systems A\ and A%, it is necessary to apply a global operation 
on the joint system AiA?. A simple global operation consists in letting systems A\ and A2 interact. In general, as a 
result of direct interactions, systems A\ and A2 become quantum correlated. On the other hand, entanglement can 
also be transferred from a third system B1B2 to AiA 2 . If systems B\ and B2 are entangled, then one can apply local 
operations on the pairs of systems (A\Bi) and (^2-62)- As a result of th ese operations it is possible to transfer the 
entanglement, originally in BiB 2 , to the joint system AiA 2 (see Fig 1(a) ). This approach is useful in the case where 
A\ and A2 represent two distant systems or when they interact weakly (or do not interact at all). The entanglement 
transfer from flying qubits to localized qubits has been extensively studied (see J3] and references therein). Entangle- 
ment transfer from two qubit systems to two qubit systems was investigated in [3|. Moreover, in Fig 1(a) the systems 
Ai, Bi, (i = 1, 2) may represent different degrees of a freedom of single particle, such as spin and momentum. In fact, 
in it was shown that the entanglement between the momenta of two different particles can be transferred to the 
spins of the particles under Lorentz transformations. In the last few years, entanglement transfer from many body 
systems and relativistic quantum field has been considered. In p|, a scheme was proposed to extract entanglement 
from a quantum gas to a pair of qubits via local interactions. Also, it was shown in [y] that two qubits interacting 
locally with a quantum field can become entangled even when the qubits remain in causally disconnected regions 
throughout the whole interaction process. The previously described scenarios are depicted in Fig 1(b) where systems 
Ai and A 2 are coupled to a common system B. 
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(a)Entanglement Transfer from 2-qubit 
system B1B2 to 2-qubit system A1A2 



(b)Generalized entanglement transfer. Here 
the operators coupling A\ to B and A2 to B 
commute 



FIG. 1: Entanglement transfer schemes. 



In the present paper we investigate the entanglement transfer scheme illustrated in Fig 1(b) where A\ and A 



represent a pair of two level systems (qubits). Here, we assume that the group of operators involved in the interaction 
between A\ and B commute with the operators coupling A 2 to B. By doing so, we mimic the entanglement transfer 
scheme described in Fig 1(a) The difference between both cases lies in the fact that the Hilbert space Hb of system 
B may not have the structure Ti.Bi ®^s 2 - This point requires further explanation; for example, when systems A\, A 2 
are coupled to different segments of a chain of coupled harmonic oscillators, then one can identify the Hilbert spaces 
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7Ybi an d Hb 2 . However, when A\ and ^4 2 couple to a quantum field vial local field operators then it is not clear what 
should be taken as the subsystems By and B 2 . On the other hand, since operators acting on different Hilbert spaces 
commute, it is clear that the situation portrayed in Fig 1(a) is a particular case of the scheme we investigate. 
The paper is organized as follows. In section (|TT|) we present some necessary conditions that the operators coupling 
the systems must satisfy in order to allow entanglement transfer from B to A\Ai- Also, we derive an expression for 
the two qubit reduced density matrix. In section (JTTTJ) , we introduce a special class of Hamiltonians describing the 
interaction between qubits and bosonic systems. For this class of Hamiltonians (being a generalization of the Jaynes- 
Cummings Hamiltonian), the reduced density matrix of the qubits assumes a particulary simple form. In section (|IVp 
we expand the entanglement measure (negativity) in terms of the coupling strength. We express the first nonvanishing 
contribution to this expansion in terms of 2-point and 4-point correlation functions involving the operators acting 
on system B. We study the entanglement of qubits in the weak coupling approximation (or equivalently, for short 
interaction times), for different N-particle excitations of the bosonic system. We also investigate the entanglement of 
the qubits in the case where system B is in a mixed state of the form ps = J2nP n 1-^0 Finally, in sections (jVJ 
and (jVIjl we compute the exact density matrix for the qubits and discuss their entanglement for different N-particle 
states and operators. 



II. FORMULATION OF THE PROBLEM 



Consider the interaction between system B, with Hilbert space H.b , and system A1A2, with Hilbert space H.a 1 ®T~La 2 ■ 
This interaction induces the following global operation on the system BA1A2: 

p^UpV^. (1) 

Assume that initially system A is in a separable state pa — ^2 n PnPAi,n ® PA 2 .n while system B is in the state ps- In 
addition, we assume that systems A and B are uncorrelated i.e. p = pa® Pb- After the interaction, the state of A is 
described by the reduced density matrix 

p A = Tr B (Vp A ®PB^) (2) 

obtained by tracing out the degrees of freedom corresponding to system B. Making use of the spectral decomposition 
Pb = YIi^b.i \^b,i) (^b,i\ and choosing an orthonormal basis {I^Bj)} for H.b, one can write the operator-sum 
representation Q of the operation pa — > p A 

P A = Y, A jIPaA^. (3) 

3,1 

The above expression shows that the positivity of the density matrix pa is preserved by the operation pa — * p A ■ The 
operators Aj.i are given by 

Aj,l = V^<*Bj|U|A fllI ) (4) 

and they satisfy the relation ( AjjAjj' = Ia which guarantees that Tr(p A ) = 1. Let U = e~ lUt where H is a 
Hermitian operator of the form 

H = 5^Ai,fc®Bi, fc +53Aa,i®Ba,i. (5) 

k I 

Here, and Bj /. are Hermitian operators acting on T~Lah * = (!) 2) and TLb, respectively. In addition, we assume 
that [B^fcjBa.i] = which implies that the evolution operator factorizes as U = Ui • U2 where U,; = e~' lt Ai < k ® Ki . k , 
i = (1, 2). It turns out that if all the operators coupled to one of the systems, say Ai, commute i.e. 



»i.fc,»i.fc' 



= (6) 



then the state p A will remain separable. In order to prove this fact, we choose {|$j)} in expression ([4]) to be a basis in 
which all the operators Bi^ are diagonal, that is, Bi.fc \$B,j) = &i,fcj \&j)- Then, the operators Aj t i assume the form 

Aj.i = Uij ® A 2 ,j,i (7) 
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with Uij = e ''^t ll ''' Au and A2,j,i = \f^B,i {® b ,j\^i\^B ,l) ■ It is clear that operators of the form (JT]) map a 
separable state into another separable state. For example, if pa — PAx ® Pa 2 then using ([J} one obtains 

PA\ ® Pa 2 ^TtMj.iPa^^jpa^ A ^PmMJ (8) 

TrA 2 (A2,j,iPA 2 A2,jy) 

which is a convex sum of density matrices i.e. p A = ~^2j iPj.iPA 1 ,ji ® pA 2 .ji, with ~^2jiPj,i = 1- Therefore, in order 
to entangle systems A\ and Ai each group of operators {Bi^} and {B^} in must contain at least one pair of 
noncommuting operators. In particular, this fact rules out operations of the form Ui = e~ l * A *® B,i . Notice, however, 
that this statement holds true for time independent Hamiltonians of the form ([5]). If one incorporates the free evolution 
of systems A and B, then, in general, the time evolution operator U(i) will contain noncommuting operators. In what 
follows, we will neglect the free evolution of the systems A and B. Furthermore, we will consider the situation where 
Ai and A 2 are two-level systems (qubits) whereas B is a bosonic system. 

In principle, one can obtain the reduced density matrix p from Kraus representation^. However, for our purposes 
it is convenient to write equation © as 

(o| / \a') = Tt b (pb (^a\V^ \a') (a\ U |^ A » (9) 

which explicitly shows that the matrix elements of p A are given by expectation values of operators acting on system 
B. Since the unitary operator U factorizes i.e. U = Ui • U2, one can define the operators = (0|U;|0) and N; = 
(1|U,|0) for (i=l,2). Assuming that the qubits are initially in the separable state <f>A = |0,0) and choosing the basis 
{|ai) = |0, 0) , |a 2 ) = |0, 1) , |a 3 ) = |1, 0) , |a4) = |1, 1)} for Ha = Hax ®7~La 21 one can write the reduced density matrix 
© as 



(Ki t KiN 2 t K 2 ) 



^ 2 f ]K 2 ) 



(Ki t KiK 2 t N 2 ) (Ki t KiN 2 t N 2 ) (N^KiKa^a) 
(Ki t NiK 2 t K 2 ) (Ki^NiN 2 ^K 2 ) (Ni^NiK 2 ^K 2 ) 
V (K^NiK^Na) (K 1 t NiN 2 t N 2 ) (Ni t NiK 2 t N 2 ) 



\ 



(Ni t KiN 2 t IK 2 

(NitNiNatKa 
(Ni^NiN 2 t N 2 ) J 



(10) 



where we used the notation 



TrsO 
tr.t 



+ nM = 



Notice that the operators 

(i = l,2) =*Tr(/) = l. 



[jtj. satisfy the relations 



(11) 



Thus, from the matrix (p~0|) we see that all the properties of state p A (in particular, separability) depend on the 
interplay of the different correlation functions of the operators {Nj, Nj , Kj, K^}. 



III. THE INTERACTION MODEL 



The most general form of the interaction between qubit Aj and the bosonic system is given by the Hamiltonian 
Hj = Ij £g> Bo,* + cr x ,i ® ^x,i + &y,i ® ®j/,t + &z,i ® ^>z.i where Mk,i are Hermitian operators. In this paper we restrict our 
discussion to the case where Bo.j = B Zi j — 0, i — (1, 2). Under this assumption, the Hamiltonian of the system can be 
written as H = Hi + H 2 where 

The evolution operator factorizes (i.e. U = Ui ■ U 2 ) if the set of operators coupled to A\ commute with those coupled 
to A2- Consequently, we impose the following conditions: 

[F 1; F 2 ] = [Fi,F 2 t] = pFit.Fa*] = 0. (13) 

Now, one can express the operators Kj = (0|Uj|0) and Ni = (1|U, |0) in terms of Fj and F^. From (fT2|) . one obtains 

Kj = (0| U, |0) = V 7^ r i 2fe (F l t F I ) fe = cos(v/F i t F i t) (14) 
f ( 2fe ) ! 



N, = <1|U< |0) = { - i )f2-L^t^W i (FM) k = -iF ^y 5 ^ . (15) 
^(2^ + 1)! VfTfT 
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Since B is a bosonic system, one may write the operators Fj in their second quantization form. We will restrict our 
discussion to the class of operators expressible as 

F(a,a^) =p{a^a)a n +q{a^a)a^ m (16) 

where p(a^a) and q(a^a) are functions of a^a and n, m > 1. Notice that while the above class includes field operators 
it does not include 1-body operators. If the eigenstates of the density matrix pb are also eigenstates of the particle 
number operator N = J2k a ^ (4'B,k)a(4>B,k) then from (fT0| . (fT4|) and (fT5|) we see that some of the matrix elements of 
p A vanish. In fact, under these assumptions p A takes the form 



P A 



/ pu pu \ 

P22 P23 

P23 P33 

\P*14 PU I 



(17) 

The eigenvalues of p A are given by Ai,± = \{p\\ + pa ± \J (pn + Pa)" 1 - 4(piip 44 - |pi4p) and A 2 ,± = \{p2 



>22 



P33 ± \J (P22 + P33) 2 ~ 4(p 2 2P33 - |p23p)- Making use of Schwarz inequality (| (x\y) \ < \\x\ \ ■ \ \y\\), one easily proves 
the relations pnp 4 4 > | P14 1 2 and P22P33 > IP23I 2 which ensure the positivity of the density matrix pa- On the other- 
hand, the partial transpose [8( of p A is defined as 

{ai,a 2 \p TAl \a'i,a' 2 ) = (a' 1 ,a 2 \p A \ai > a 2 ) 

and its matrix representation is 

(pn p* 23 \ 

T Al = P 22 P\, (lg) 
I P14 P33 I 
\P23 P44/ 

Notice that p rA i can be obtained from p replacing P23 by p^ 4 and P14 by P23. Hence the matrix p T/1 i will have a 
negative eigenvalue if either 

"-23 = 1/023 1 2 ~ P11P44 > or "14 = |pi 4 | 2 - P22P33 > 0. (19) 

If one of the above inequalities is fulfilled, then according to the PPT criterion [§], the two qubit system ^1^2 will be 
in a nonseparable (entangled) state. Furthermore, in this particular case, one can show that p A > implies that 7123 
and ni4 cannot be positive at the same time. For example, suppose that 77,23 > 0, then one can write the inequalities 

P22P33 > \P23? > PnP44 > Ipul 2 

which imply 7144 < 0. Therefore the partial transpose p TA i can only have one negative eigenvalue. In fact, it can be 
proved that in 2 x 2 dimensions, the partial transpose of any density matrix can have at most one negative eigenvalue 

In the previous section, we assumed that the initial state of system A1A2 was \4>a) — = 1 0, 0) . Consequently, 
we wrote down the matrix (|10p representing the final state of the qubits. However, if system A1A2 is initially in 
a different state then, obviously, the matrix elements of p A will be different from those in expression (fT0| . Let us 
consider, for example, the case \4>a) = |1, 1) = f x (8 cr x 1 0, 0) . Now, equation ((9]) can be written as 

(a| p A \a') = Tr B (p B (ai| o*JJ[a x ® \a') (a| 4^1^ ® alU 2 (r x |ai)) (20) 

where |a) is obtained from |a) by flipping both qubits, i.e. |a) = a x <8> <J X \a). Since 



riUka. = (^t »j, (21) 

it is easy to conclude that if initially system A was known to be in the state {0,4} = |1, 1) S TLax <8 then the final 
density matrix p" 4 is given by 



'p 44 p 4 



V pi4 pn, 
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where the matrix elements py are determined by means of the expressions (jT0|) , ([14]) , (Tl5|) with the operators replaced 
by the operators F^ and vice-versa. Treating p A as a function of the initial state 10^) and the operator U, one can 
summarize the previous discussion as 

p A {\a A ) , U(F fe )) = V P A (\ ai ) , U(F+))W (23) 

where V = a x ® a x = \a\) (a<i| + {0,4} (a\\ + \a 2 ) (^3) + | CI3} (02) is the unitary transformation corresponding to the 
basis permutation (1, 2, 3, 4) — » (4, 3, 2, 1). Similarly, one can write 

p A (\a 2 ) ,U(Fi,F a )) = V/flai) ,U(Fi,f£))V+ (24) 

with V = / ® <r x . 

IV. SERIES EXPANSION FOR JV(p A ) 

In general, for a given set of operators {Fi,F2j-, it may not be possible to determine the matrix elements of p A in 
closed form. On the other hand, it is interesting to explore the different terms in the series expansion of 7123, n i4 ( see 
IT9"| in terms of the operators {Fj,Fj}. From the general expression (TIT)]) for the reduced density matrix, we expect 
this series expansion to contain correlation functions of the operators {Fj,F2^} . Introducing a coupling constant, i.e. 
Fj — * g{¥i, one obtains a series expansion of the form 

n(g\,g 2 ) = niglgl + n 2 (gfg% + g\g\) + ... (25) 

where n denotes either n 2 3 or 7114. The above form can be justified as follows. Each term of the series expansion 
should be symmetric in gi and g 2 . Terms containing odd powers of gi, e.g. gig 2 + 3i<?2, should not be present since 
then one could change the sign of n by simply changing the signs of g\ or g 2 . Clearly the entanglement in A\A 2 
should not depend on the sign of g\ or g 2 . Finally, if we switch off one of the interactions (gi=0 or g 2 = 0) then both 
rt23 and nu vanish (see expressions (|19p). This rules out terms of the form gf + g 2 . Now, we proceed to express the 
first nonvanishing contribution to the expansion (|25p in terms of Fj and Fj'. Making use of the series expansion (|14p 
and Q3) for K 4 and N h one writes K t = I - ^F^F*)* 2 + ^(F^Fj) 2 * 4 ... and = -i¥it(I - -^F^F^ 2 ) + . . .. From 
(I10p one obtains the following series expansion for the diagonal matrix elements p^: 

pn=I- (Pi) t 2 - (P 2 ) t 2 + ~«P 2 ) t 4 + (F 2 ))t 4 + (P X P 2 ) t\ p 22 = (P 2 ) t 2 - i (P 2 ) t 4 - (PiP 2 ) i 4 , (26) 
P33 = (Pi) t 2 - i (P 2 ) t 4 - (PiP 2 ) t 4 , P44 = (PiP 2 ) t 4 (27) 

where P., = (F^F,). Similarly, one expands the P23 and P14 in series of {Fj,Fj^} obtaining the following expressions: 
P23 = (Fi t F 2 )t 2 -i(F 1 t F 2 t F 2 2 )t 4 -i(Fi t2 FiF 2 )< 4 -i(F 1 t FiFi t F 2 )< 4 

_i.( Fl tF 2 F 2 tF 2 )f 4 (28) 
(Fi t F 2 t ) t 2 + \ (Fi t F 2 t2 F 2 ) t 4 + \ (Fi^FiFa 1 ) t 4 + 1 (Fi t F 2 t F 2 F 2 t ) t 4 



Pi4 = -(Fi T F 



3! 



+ (F 1 t FiFi t F 2 t )i 4 . (29) 



Thus up to fourth order one has 



Pupa = (P1P2) t 4 = (Fi t FiF 2 t F 2 ) t 4 = ((F 1 F 2 )tF 1 F 2 ) t 4 > (30) 

|pi 4 | 2 = I (Fi t F 2 t ) \ 2 t 4 (31) 

P22P33 = (Pi) (P2) t 4 = (F^Fi) (F 2 +F 2 ) i 4 > (32) 

|p 23 | 2 = I (Fi f F 2 ) | 2 i 4 . (33) 
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For consistency, one can check the density matrix is positive. In fact, the identity (AA^) > | (A) | 2 with A = Fi^F 2 ^ 
implies that pupa > \pi4\ 2 while the inequality P22P33 > I/O23 1 2 follows from Schwarz inequality. Clearly, the partial 
transpose p TA i may be negative. In fact, substituting the above approximations for p A in the expressions (fT9j) for n 2 3 
and ni4 one obtains: 

n 23 = t 4 (\ (Fi+F a ) | 2 - (F 1 tF 1 F 2 tF 2 » (34) 
n u = t\\ (F 1 t F 2 t ) | 2 - (Fi+Fi) (F 2 t F 2 )). (35) 

For small values of t, the above expressions can be used to detect the presence of entanglement in the system A1A2. 
If one of the above quantities (n 2 3 or nu) is positive then the state of qubits A\ and A2 is nonseparable. Notice that 
if we choose the operators F^ to be field operators acting on B, then n 2 3 and nu will contain 2-point and 4-point 
correlation functions of these operators. On the other hand, if the operators {Fi,F 2 } are normal (i.e. [F^F^] = 0) 
then n 2 3 and nu are be negative. This is a reflection of the fact that in order to entangle system A\ with system 
A2 we must have non commuting operators. In the particular case where the operators F^ are linear combinations of 
creation and annihilation operators we have 

[Wi,¥l]=Ci for (* = 1,2), (36) 

where Cj is a c-number. In this case one can show that n 2 3 and nu are negative if either c\ < or c 2 < 0. For 
example, let c 2 < 0, then Schwarz inequality leads to 

nu = i 4 (| (FiW) I 2 - (Fi^Fi) (F^) + c 2 (F^)) < c 2 t 4 (F^) < (37) 
n 23 = t\\ (F x tF 2 ) | 2 - <(F 1 t F 2 )(F 1 t F 2 ) t ) + c 2 (F^Fi)) < c 2 i 4 (F^Fi) < 0. (38) 



The relation (|38p is obtained using the inequality 

(AA+) > I (A) | 2 . (39) 

Finally, the entanglement measure we use in this paper is the negativity N(p A ) defined as twice the absolute value 
of the negative eigenvalue of p TA i [ll|. In our case, the eigenvalues of the partial transpose p TA i (fT8|) are (up to fourth 
order in t): 



V(Pn + Pa) 2 + 4 ™2 3 ) « r — w -n 23 (40) 

Pll + P44 

V(P22 + P33) 2 + 4n 14 ) « ^— . (41) 

P22 + P33 



\ Ta i l ( 



P44 



Xo Al ' 



^(P22 + P33 



A. Examples 

Utilizing expressions J34|) and (|35|) one can calculate the entanglement transferred from system B to the pair 
of qubits AiA 2 . In this subsection we compute the quantities n 2 3 and nu for interactions of the form (|16p i.e. 
Fi(a,a)) — f(a)a)a n + g(a>a)a) and statistical mixtures of N-particle states pb = ^2nPn \N) (N\. Most of the 
examples that we present here can be computed exactly using (|10j) . Nevertheless we believe that the first nonvanishing 
contribution to entanglement gives us some hints about the states and interactions that induce entanglement in system 
A\A2- Moreover, from a technical point of view, expressions of the form (F]/F 2 ^), (F;jFiF 2 ^F 2 ) . . . can be easily 
computed using Wick Theorem [l5[ . 

First, we consider the JC [l2[ describing a system of two two-level atoms interacting with a set of electromagnetic 
modes. Let the creation and annihilation operators for these modes be and at. We assume that the Hamiltonian 
for the two atoms is of the form HI = Hi + H 2 with Hj, (i = 1, 2) given by 

0,\ [0 1 



= ^2gt,k(<7+a k +a-a k r ), er+ = f 1 Q ' j , = L A . (42) 
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In this representation, |0) = denotes the ground state of the atom while |1) = ^\ represents the excited state 
of the atom. Let \k) denote the mode created by the operator at. Introducing the state |0j) = J\ 7=^ , \k) we 



rewrite the Hamiltonian (l42l) as 



gi(a + a(4>i) + o-a)(4>i)) = g t ( 



ot(0i) 
a(<t>i) 



(43) 



where <?.; = ^J^2 k gf k - This Hamiltonian is of the form H, = ( ^ ^ J with Fj = gi<d{<j)i). These operators satisfy 

[Fi,F 2 ] = 0, [Fj,Fj] = gf and [Fi,F 2 ^] = 3152 (01 1 ^2) • If we assume that the states \<px) and |0 2 ) are orthogonal, 
then relations (fT3|) will be satisfied and we can apply, the results from previous sections. From the physical point 
of view, the orthogonality of \4>i) and | ^2) could describe the situation in which two atoms are sensible to different 
modes i.e. they absorb/emit particles in different modes. If initially the qubits are prepared in the state \4>a) = |0, 0) 
while system B is in the state = -^(a^ (4> B )) N |0) (representing N bosons occupying the same state |0b)) one 
easily finds that 

(FitFa) = 3192 (ai f a 2 ) - 9W2N (0 s |0i) (falfo) (44) 
(Fi t FiF 2 t F 2 ) = ghl{aSaWa 2 )=gtglN(N-l)\ (^ B ) | 2 | (0 2 |0 B ) | 2 . (45) 

Substituting these results in expressions (f34|) and ([35]) one obtains: 

^23 = Ng 2 g 2 t 4 \ (falfa) | 2 | (0 2 |0 B ) | 2 > (46) 



n 



U = -iV 2 2 2 * 4 |<0i|0 B )| 2 |<0 2 |0 B )| 2 <O 



From the above expressions we conclude that system A\A 2 is always entangled except for the cases where u\ = 
(01 1 0s) = or u 2 = (02 1 0b) = 0- One can also consider the situation where the initial state of system B is a mixed 
state of the form 

00 00 
Pb = ^2pn\N)(N\, X>" = 1 ' (47) 

iV=0 N=0 

and \N) — "y^y( a ^(0s)) Ar |0). As discussed in section all the expressions we have derived so far hold for the class 
of density matrices having eigenstates with fixed particle number. Thus, from equations ((44)) and (|45|) . one obtains 

00 00 

^23 = S 2 9 2 2 * 4 M 2 M 2 ((£ P^) 2 - £ P^iV(iV - 1)) (48) 

JV=0 iV=0 

which implies that p A will be nonseparable for probability distributions pm satisfying 

00 00 

(J2pnN) 2 >J2pnN{N-1). (49) 

N N 

It is interesting to_note that the above condition implies that the probability distribution {pn} must be sub-poissonian 
[Tij .i.e. N > N 2 — (N) 2 = er 2 ^. An example of a distribution satisfying TV > ajj is the binomial distribution 

Pn = {n)p N ( 1 ~ P) M ~ N ( for N = 0, 1, ■ ..M) having N = Up and a 2 N = Mp(l - p). In this particular case one 
obtains 

"23 = g 2 g 2 2 t 4 K| 2 M 2 A/p 2 > 0. (50) 

Finally, notice that a thermal-like state p B = (1 — 2) ^2 N zN 1-^0 (-^1 yields a separable state for A\A 2 . If the initial 
state of the two-level systems is \4>a) = 1 1 , 1) then according to section . the entanglement of A\A 2 is determined 
by n 2 3 and 7114 with the operators Fj being replaced by F^. As expected, in this case we obtain a separable state. In 
fact, if we denote u, = (0i|0 B ), then we have 

I (FiFat) | 2 - (FiFitFaFat) = -g 2 g 2 N(\ Ul \ 2 + \u 2 \ 2 - \ Ul \ 2 \u 2 \ 2 + 1) < 0. (51) 
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One can also consider a state of the form 

= VJVi 1 , jV2 , ( at (^)) JVl (a t (^J) JV2 |0> = \m,N 2 ) (52) 

describing the situation where N± bosons occupy the state \4>b x ) and N 2 bosons occupy the state \4>b 2 ) orthogonal to 
\4>Bi), i-e. {4>b 1 \4>b 2 ) — 0- Straightforwardly, one finds that 

Oi|$n) = V^ , «4,i|JVi-l,JV 2 ) + -v/^'tii ) 2|JVi,JV2-l) (53) 
I (F 1 t F 2 ) | 2 = 5 2 52 2 |^ lU * iaU24 + JV2U* ll2 U2,2| 2 (54) 



where tij,fc = (<pi\4>B k )- On the other hand, (Fi^FiF 2 ^F 2 ) reads: 

(Fi+FtFa+Fa) = ff?s§(Wi(Wi ~ l)|u 1)lU2 ,i| 2 + N 2 {N 2 - l)\u 2 . 2 u h2 \ 2 + A^Ki^o + u 2 , 2 u 1A \ 2 ). (55) 

Combining the above results one obtains 

"23 = glg 2 t 4 {Ni\ui,iu 2 ,i\ 2 + iV 2 |wi,2U2,2| 2 - NiN 2 (\u 2 ,iui, 2 \ 2 + |u 2 , 2 wi,i| 2 )). (56) 

Here, two things are worth mentioning. First of all, equation (|56[) indicates, that in this approximation, no en- 
tanglement will be transferred to A\A 2 for large values of N\ and N 2 . Also notice that when Ni = N 2 = 1, 
| (Fi^F 2 ) | 2 = (? 2 (? 2 |iViu*i.iu 2 ,i + N 2 u* i j2 w 2j2 | 2 vanishes when the 2x2 matrix Ui r k = (4>i\4>B k ) is unitary implying 
n 23 < 0. This case will be studied in detail in section (|VI[) . 



B. Algebraic Construction of Operators 

One can also consider the entanglement induced in system A\A 2 when the operators Fj are linear combinations of 
creation and annihilation operators. Let 

Vi = g i (ai + \p i \e iei a^), i = (1,2), with [01,03*] =0 and < 1. (57) 

If system B is in a 1-particle state a B |0), then 

{¥!%) = (Q\a B Vi%a B i\0) = gi 9 2{u l *u 2 + Ul u 2 *\0 1 \l3 2 \e- i ^- 8 ^) (58) 
(Fi t FiF 2 t Fa> = (0|a B Fi t FiF2 t F 2 aB t |0) 

= glg 2 2 (\Pi\W + Kl 2 (i + l/3i| 2 )l/3 2 | 2 + M 2 (i + lft| 2 )|/3i| 2 ) (59) 

In the case where u\ — u 2 — -j= and (3\ = f3 2 = j3 we have 



(Fi+Fa) - (Fi+FiFa^a) > for \0\ < \/3 max \ = J 2 ^ 1 = 0.51 (60) 

and n 14 = ^((F^F^) — (F^Fi) (F 2 *F 2 )) < for f3 € (0, 1). Thus, one can mix creation and annihilation operators as 
in equation (|57)) . obtaining an entangled state for A\A 2 for < \(3 max \. 

We conclude this section by constructing a set of operators inducing entanglement in system A\A 2 when B is the 
particle vacuum state |0). This problem was studied in [6j where it was shown using perturbation theory that a two- 
level systems can become entangled after having locally interacted with a scalar field. Here, we assume an effective 
interaction between the qubits and system B of the form 

H 4 = (S. F (Q with ¥, = g^afa) + Ptat^i)). (61) 

We need to accommodate conditions (fT3"|) into this picture. Notice that if {4>i\(f> 2 ) = (ipx\ip2) = then [Fi,F2*] = 0. 
On the other hand, [Fi, V 2 ] = holds if we impose the condition 

ft (V>i|<fe> = /MWi) • (62) 
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Now, we have 

(iV|F X tF 2 t|iV) - 5l52(/3l(V'l|02)+A r ((0B|0l)/?2(V'2|0B> + (0B|02)/?l(V'l|0B» (63) 

(7V|FM|iV) = ,g 2 2 (|A| 2 + N(\ (4>b\^) I 2 + m 2 \ (4>b\^) I 2 ))- (64) 

If \§b) = |0), we set N — in the above equations and compute nu vac from (|35[) obtaining the expression 

ni4vac = 5 2 5 2 2 i 4 |/3i| 2 (l (^i|02> | 2 - |/? 2 | 2 ) (65) 

which indicates that the system A\A 2 is entangled for < |^| < | ("01 i 02) I ( or equivalently, < \(3\ \ < \ (ip2\4>i) |)- 
We close this section by presenting a situation in which operators of the form Fj = gi(a((j>i) + /3*(r(ipi)) appear in the 
effective interaction between the qubits and the bosonic system. Consider two qubits interacting with a relativistic 
scalar field [l5| . Let the Hamiltonian of the system be 

M = H + J2 9i(t)<Tx,i4>i, <Jx = (J J) (66) 

i=l,2 ^ ' 

where Ho = \ 2i=i 2 ^A^z.i + 53fe w fc( a fc^ a fe + 5) i s the free Hamiltonian of the qubits+field system. The functions 
<7l(t) and 52 (i) describe time dependent coupling strengths and 4>i — f Q d 3 rfi(f)4>(t,r), <j> 2 = J Q d?rf2{r)4>(t,r) are 
average fields on the spatial regions 0\ and O2 where the qubits are located. If the interaction is fast compared to the 
free evolution of the qubits and the field (by taking the average fields we introduce a cut-off for the field frequencies) 
one can make use of Magnus approximation jT(| and write the time evolution operator as U = e -i $ dm,(t) w h ere 
H/(£) is the interaction picture Hamiltonian. Following, [f| we assume that the qubits remain causally disconnected 
throughout the whole interaction process. The time evolution operator can be written as U = e~ lffc// where 

H »ff = E / (A« e T") = E {I F o ) • (67) 

Here the operators Fj read 

Wi= / dtgitye^Mt) = I „ * h{k){gi{uk - ^)e^a fc + g*(w k + u Ai )e~ ik - ri a\) (68) 

JO J (2n)2 y/2LJ k 

where §i{uS) is the Fourier transform of the function gi(t) i.e. g"i(u)) — J dte l( ^ u ' t gi{t) while fi(k) is the Fourier transform 
of the smearing function ft (r) . For a massless field and qubits with u> Al — lua 2 , one can prove that in the limit where 
the volumes of the regions 0\ and O2 approach zero, one has 

WiM}= I (n $^= {\~9i^k-^ Ai )\ 2 - \§i(Uk +"A i )\ 2 ) = ^ [dtg*(t)>0, t = (1,2). (69) 
J (27r)2^2w fe 27T J 

Thus the necessary conditions (|36p for the operators Fj are satisfied. Notice that if we neglect the free evolution of 
the qubits (we set ujai = ) Fi and ¥2 become normal operators and the state of the qubits remains separable after 
the interaction. The final entanglement of the two qubit system depends on the form of the functions <?«(£), i = (1,2) 
and it will not be discussed here. See \m for a detailed discussion. 



V. ENTANGLEMENT FROM N PARTICLES OCCUPYING A SINGLE 1-PARTICLE STATE. 

In the previous section we found that for interactions of the form (TT2"|) with F, = gia(<f>i), the system A1A2, initially 
in the state, |aj.) = |0,0), becomes instantaneously entangled when B contains N particles occupying a single 1- 
particle state \4>b)- The only requirement for the state \<j>g) is to overlap with \(j>i) and |^ 2 ) i- e - ui = (0i|0s) 7^ 
and U2 = (</>2|0b) 7^ 0. Since the states \<j>\) and \<j>2) are orthogonal, any 1-particle state can be written as \<Pb) = 
J2i u i a \ |0) + ut \4>t) where (&|0r) = for i = (1, 2). Therefore the state |$g) = y^n ^- 8 )^ 1^) can ^ e written in 
occupation number representation as 

|*B> = £ ]l nM (N N -n 1 -^ 1 ^- ni - na \—- N -^- »»> <*» 
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Now making use of (fl4|) . ([15]) and (flO]) . we compute the reduced density matrix for A\A 2 as a function of time. For 
example, the diagonal element 



p n = (k|KiK 2 ^K 2 ) = <$ B | cos^^PV^t) cos 2 ( V /F 2 tF 2 t)|$ B ) (71) 

can be easily computed using the occupation number representation (|T0[) for the state Likewise, one may 

determine the remaining diagonal entries of p A which can written in the compact form 

pkk= 2. ni!n2!(JV - ni - n3)1 J> n 1 ,^(«i,m)i^(n 1 ,n a ) (72) 

"1,12=0 v ' 

where P ni , n2 ( Ul , u 2 ) = \ Ul \ 2n > \u 2 \ 2n * (1 - | Ul | 2 - \u 2 \ 2 ) N ~ n ^ and 

Fn(m,n 2 ) = cos 2 ( % /rugii) cos 2 (^/n 2 "g 2 t), F 33 (ni, n 2 ) = sin 2 ( v /nT5'ii) cos 2 (^/n 2 "g 2 <), (73) 
F 22 (rii,n 2 ) = cos 2 ( 1 /niffit)sin 2 ( v /n 2 "5 2 i), F 44 (m, ra 2 ) = sin^yroT^) sin 2 ( v /n 2 "g 2 t). (74) 



From equation (JTUJ) , it can be easily seen that operators of the form Fj = gi<ii lead to P14 = 0. Hence, the only 
non-vanishing off-diagonal matrix elements are p 23 = (<1>_b|Ni^KiK 2 ^N 2 |$b) and, obviously, p 32 = p 23 *. Again, using 
([70]) . one obtains 

^ N\ 

p 2Z =u\u 2 > — : — — ^.Q nu n 2 { u uU 2 )F 2 s(ni,n 2 ) (75) 

' m!n 2 !(iV - m - n 2 - 1)! 

where the polynomial Q ni ,n 2 (ui,u 2 ) = |ui| 2ni |u 2 | 2 " 2 (l - |ui| 2 - \u 2 \ 2 ) N - ni - n *- 1 and 

jp / , , , — ,,sin(Vm + Igit) — sin(Vn 2 + lg 2 t) 

F 2 3{ni,n 2 ) =cos(V^i5i^) , , = cos{^/r^g 2 t) . (76) 

V^l + 1 V "2 + 1 

Unfortunately, it does not seem possible to express the sums (|72[) and ([75]) in closed form. However, numerical analysis 
indicates that one obtains entangled states for the qubits for all values of N with larger values of negativity when N 



is an odd number (see Fig 2(b) ). The above results can be extended to the case where 

¥i= gi af, t=(l,2), m>l, and [01,02*1=0. (77) 

This type of operators describe the situation where to particles are needed to flip (or excite one of the atoms) the 
state of one of the qubits [l7|. The expressions ([72"]) and ([75]) (corresponding to the case m=l) can be generalized 

to any value of m. Using the identities ata^™ — ^7^T~yp^ > o-i m o-i — Jt' a .-m)! ano - writing the state as in 

equation (|70p . one arrives at the following expressions: 

ENl M 
— : — (Ttt ^P nun2 (ui,u 2 ) Fkk{ni,n 2 ), (78) 
ni!n 2 !(iv — ni — n 2 )\ 

„ N\ (m) (m) 

P23(™) = K«2) m > — : — Vflvr tt Q ni .„ 2 (ui,u 2 ) F 23 (ni,n 2 ) (79) 

' n\\n 2 \[N — ni — n 2 — m)\ 

ni,ri2 

(m) (m) 

where F kk (n u n 2 ) = F kk (j^L^, j^^) , Q „ 1>n >i,u 2 ) = H^H 2 "^ - Kl 2 - |h 2 | 2 )^-^-™ and 



(«>, , «»( Vt^^') sin (\/ i! ^5i0 cos( v /3I 52 i) sin( v /I-^ 52 

^23 (ni,n 2 ) = ? , ? , ■ (80) 

[nx+my. / (n 2 +m)! 



«2 ! 



From equation (|78p one finds that p 44 vanishes for states having N = {to, to + 1, ... , 2m — 1} particles occupying 
the 1-particle state \4>b)- This is easy to understand; the matrix element p 44 corresponds to the process in which the 
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Negativity, N=Odd 
1.0 r 




Negativity, N=Even 
1.0 r 



0.8 
0.6 h 
0.4 
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(a)N=l (dashed line) and N=3 (solid line) 



0.2 



,A\ ,„,/ A, 



' v., !\ 



2 4 6 8 10 

(b)N=2 (dashed line) and N=4 (solid line) 



FIG. 2: Negativity as a function of time for odd and even values of N. We have assumed gi — g2 — § and ui = ui = 



states of both qubits are flipped. Therefore, the presence of at least 2m particles is required to have P44 7^ 0. In our 
model, /544(f) = implies that the two qubits are entangled with negativity N{p A ) = \J {pii) 2 + 4|p23 1 2 — P11 > 
for any value of t (except for those values of t for which /023(f) = 0). For N = m, \ui\ = 1 1*2 1 = ^ and gi = <?2, 
expressions ([78")) and ([75)) take the particulary simple form 



P23 = sin 2 (y/W.gt), p k k = ^ 

n\+n2=N 



(m) 

F kk(ni,n 2 ) 

m\n 2 \ 



for fc=(l,2,3) and /9 44 = 0. 



(81) 



Here, we notice that the maximum values of the negativity Af(p A ), behave like 1/2 2N 1 for large values of N. For 
m > 1, this behavior can be improved. In fact, for TV = m + 1 we have p 2 % = ^Mi sin(Vm\gt) sm(y/ (m + \)\gt), 
P44 — and as a result Af(p A ) ~ N/2^ 2N ~ 3 ^ for large values of N. Graphs for the cases m = 2, N = 2 and m = 2, 
N = 3 are shown in Fig 3(a) In what follows, we will restrict our discussion to the case m = 1. 



A. Entanglement from 1-Particle States 

If the qubit system is initially in the state |0, 0), we obtain from ([72")) and ([75)1 the following density matrix 
p n = l- |ui| 2 sin 2 (5it) - \u 2 \ 2 sin 2 (g 2 t) p 2 2 = \u 2 \ 2 sin 2 (g 2 t) 



P33 = K| 2 sin 2 (gi<) 



P23 = ui*u 2 sm{gxt) sin(g 2 t). 



As we know from previous discussions, this state is always entangled with negativity J\f(p A ) — \J (pii) 2 + 4|p23 1 2 — P11 ■ 

i and gi — g 2 , the negativity oscillates between zero and 



"2 



Furthermore, if the amplitudes Ui satisfy 

its maximum value Af = 1 with period — (see Fig 2(a) |. However, this behavior changes substantially if the initial 



state of A\A 2 is \a±) = |1, 1). From section . we know that the reduced density matrix can be obtained from the 
relation 



/(|a 4 ),U) 



(P44 p 41 s 

P~33 P32 

P*32 P22 

\pl4 pn, 



Vp A (\a 1 ),U(¥ 1 \¥ 2 ^)V\ V 



(82) 
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In this case we also have pu = and therefore in order to quantify the entanglement in system A\A 2 one needs the 
following matrix elements: 

p 4 4 = cos 2 (git) cos 2 {g 2 t) + \ui\ 2 cos 2 (g 2 t) (cos 2 (V2git) -cos 2 (g 1 t)) (83) 

+ \u 2 \ 2 cos 2 ( 9l t)(cos 2 (V2g 2 t) - cos 2 (. 92 i)) 

Pll = sin 2 ( 9l t) sin 2 (g 2 t)-\ Ul \ 2 sin 2 (g 2 t) (cos 2 (V2 gi t)- cos 2 ( 9l t)) (84) 

- \u 2 \ 2 sin 2 (g 1 t)(cos 2 (V2g 2 t) - cos 2 (g 2 t)) 

Pzi — u\u\ cos(V%g 2 t) cos(y / 2ffii) sin(git) sin(g 2 t). (85) 



The above expressions dictate the time dependence of the negativity (see Fig |3(b)[ ) and contrary to the situation in 
which the initial state of A\A 2 was |0, 0), now system A\A 2 exhibits periods of entanglement death and entanglement 
revivals. Nevertheless, for the symmetric scenario ( g\ — g 2 = g and = \u 2 \ = -^) one can obtain an almost 

maximally entangled state (AT « 1). Notice from (|83l l84l l85|) that one can simultaneously have pn « 0, P44 ~ and 



P23 



when 



V2gt : 



and gt w (2m + 1)^ => (2m + 1) w nV2. 



(86) 



The first three pairs of numbers of the form (n, 2m + 1) satisfying (approximately) these equations are (5,7), (12,17) 
and (29, 41), (5j/2 = 7.07, 12^ = 16.97, 29\/2 = 41.01). The first two pairs correspond to the peaks with 
Af « 1 in Fig|3(b)| At these points the system A\A 2 is in the state p A fts \ip + ) (ip + \ with = \ (\a 2 ) + I03)). 



Negativity m=2 

l.Or 



Negativity 

l.Or 




(a)Cases N = 2, m = 2 (solid line) and N = 3,m = 2 (dashed 
line). Here |«i| = \u 2 \ = -4= and gi = g 2 = 1. 



0.8 - 



0.6 - 




5 10 15 

(b)Negativity for \4>a) = 1 1 , 1 } - Here |ui| = |tta| = -4= and 
gi = g 2 = %. The maximum value found is M = 0.98. 



FIG. 3: Negativity vs. time. 



B. Operators of the form Fj = gt(ai + \/3i\e % *ai ) 

Another type of operators Fj belonging to the class (fTB|) are the operators of the form Fj = <7j(dj + \(3i\e l6i aJ) 
with |/3j| < 1 and [ai,ci2^] = 0. For simplicity, we assume that system A\A 2 is initially in the state |0, 0). For these 
operators, one can sum the series (fT4")l and (JTSJ) by means of a Bogoliubov transformation [l8|]. In fact, one may set 
f3i = tanh(rj) and write Fj = gibi where 6j = cosh(rj)aj + e l0i sinh(rj)aj^ and gi = cos ^ r .) ■ Using the identities 

b i =U i a i u}, with w i = e 5(^*"i 2 -W 2 ) ) z = ri e ie * (87) 

one can express the vacuum corresponding to the operators Oj in terms of the eigenstates of the operators rij = 6j 6j. 
That is 

|0) = l 2n i' 2n ^ °> ' c »* = / \, , e '"' 9 ' tanh"'(r ^ x ' 2 "'' " 



ni ,ra 2 
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In this new basis the 1-particle excitation |$s) = a) B |0) with u\ = u 2 = assumes the form 



(89) 



**> = "4 E ^^( V2 ?, + ! I2»i + l,2n 2 ) + V2 " 2 + 1 |2m,2n 2 + 1)). 
V2 ^ ' coshm) cosh r 2 

v ni,n2 

;enstates of the operators F^Fj with 



Assuming ^! = g 2 and taking into account the fact that the states |rij) are eig 
eigenvalues gfrii we obtain: 

/ AB —EF \ 

I o ^ b!±f! o 1 



o ° 



(90) 



V —EF CD J 
with A,B,C,D,E,F given by the following series 

A = Yl ^TtA l c "| 2 cos 2 (V2n+l£t), fl = V |c„| a cos a (V2n$t) 

z — ' cosh (r) — ' 

C = E 2 ",tA l c "| 2 sin 2 (V2^TT^), D = V \c n \ 2 sin 2 (V^gt) 

» cosh (r) V 



(91) 
(92) 



and 



£ = 



F = 



E 



V27 



1 



cosh(r) 
sinh(r) 



E 



rl | 2 cos(V2ngi) sin(V2n + Igt) 
2n + l 



cosh(r) 2 ^ y/2(n + l) 



|c„| 2 cos(V2n + lyi) sin(v/2(n + l)gi) 



(93) 
(94) 



Notice, that in the limit |/3| = 1, one has [F,, FJ] = 0. Hence, according to section . the entanglement in system 
A should disappear as we approach \(3\ = 1. In FigU] we present graphs of entanglement versus time for different 
values of \/3\. Numerical analysis shows, that entanglement is more strongly deteriorated for values of (3 close to 1. 
Thus we conclude that operators being mixtures of the form a + (3a* with (3 < 1 can also transfer a substantial amount 
of entanglement to the two qubit system. 



Negativity 




FIG. 4: Negativity vs. time for different values of (3. Here we assume gi = gi = 1 and u\ — U2 = The cases /3 = 0, /3 = 0.5, 
and (3 = 0.7 are represented by dashed, solid, and thick solid lines respectively. 
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C. Entanglement from Mixed States 



The expressions (f72j) and (|75|) can be also used to determine the final state of the qubits when system B is in a 

N' 



mixed state of the form ps — ^2nPn \N) (N\. In this case, the matrix elements of p A read 



4 = 5>WA0 (95) 

N 

with Pij{N) given by equations (|72l) and (jTB"]) . The separability of the qubits A\A 2 depends on the distribution {pn}- 
For example, consider the binomial distribution p^ — (^)p N (1 — p) M ~ N (which was already discussed in section PV|) ). 
In this case, equations (|95|) . ((72|) and (j75j) yield 

p** = E i „ M M! ^(^luil^^l^D^a-Cv^l"!!) 2 -^^!) 2 )^"" 1 """^^!.^)- ( 96 ) 

ni'.ri2'-(M — ni — n->)\ 

ni,n 2 v ' 

Likewise, we compute p 2 ^ from (|75p to find that the matrix elements p^- corresponding to the binomial distribution 

Pn = (x)p N (l - p) M ~ N have the same form as the p io (N) from and (J75J) with N replaced by M and the 
amplitudes u, replaced by y/pui- Therefore, this case reduces to the previously studied situation where ps = \M) (M\. 
On the other hand, since entanglement disappears as the amplitudes Ui = {(j>i\<p B ) approach zero, one expect to 

obtain a separable state p A for a Poissonian distribution i.e. ps = Pn \N) (N\ with pn = j^ts~ x (recall that Poisson 
distribution is the limit of the binomial distribution for M — ► oo, p — ► and Mp — > A). In fact, for a Poissonian 
distribution, one obtains from (19511 (l72l) and (1751) the state 



P A = 



/cic 2 \ 

c\S2 m\*m 2 

m\mi* c 2 s\ 

V Sl s 2 ) 



(97) 



where the functions Cj, s^, and m, given by 



Cj = e -A|«*l 2 V^pico S 2 (^*), Sl = e^-l 2 ^^psm 2 (VH-^), 
/r \u, i 2 (Aliii | 2 ) Tli . , . sin(\/nj + lad) 

A matrix with the structure of (f9"T| must necessarily be separable. Notice that the entanglement condition \p 2 z\ 2 > 
pupa is not compatible with the positivity condition | /O23 1 2 < p 22 P33- Hence, p A is separable. 



VI. ENTANGLEMENT FROM N PARTICLES OCCUPYING DIFFERENT 1-PARTICLE STATES 

So far we have considered N-particle excitations of system B with all the particles occupying the same 1-particle 
state. It is also interesting to study multiparticle states of the form 

N 

\<S> B ) = l[aH<j>B k )\0), {4>B k \<t>Bw)=h,v (98) 

k=l 

representing N identical particles occupying mutually orthogonal 1-particle states. Again, we assume that the interac- 
tions between B and Ai are of the form (|12[) with F< = gia(4>i) and (<f>i\(/>2) = 0. It is clear that now the entanglement 
transferred to AiA 2 depends on the relative geometry of the set of states {</>b 1 ,0b 2 , ■ ■ ■ ,4>b n } and the states |0i) 
and |02)- One can find the density matrix for states of the form (|9"8")l using expressions (fT2j) and ([75]) . Taking the 

linear combinations \4> B ) = Y,k=i x k \<pB k ), \4>b) = Y,k=iVk \4>B k ) and defining the states = ^7=7^ (<Ab) |0), 

\<&b) = -^==a^ ' (4>b) |0), one computes the auxiliary matrix element 

pn(x,y) = ($ B |K 1 t K 1 K t 2 K 2 |$ s ) = V , -P nun2 (x,y)F 11 (n 1 ,n 2 ). (99) 
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In the above expression, the polynomial P ni ,n 2 {x,y) is given by 

P ni) „ 2 (a;,y) = {mu^ {u 2 u* 2 ) n * {$ B \<t> B ) - Ul u\ - u 2 u* 2 ) N ^- n \ (100) 

with Ui = "^frXkUi^k, u* = J2kVkU* k and Ui t k = {4>i\4>B k )- From the above expression we can extract the first 
diagonal element of the two qubit density matrix. Notice, that pn (which corresponds to the original state \$>b) = 
a ^(0Bfc) |0)) is related to p~n(x,y) as follows: 

1 d 2N 

P n = M« a a a Pl,l(x,y) 

N\ ox\ . . . oxNoyi . . . oyN 

1 d 2N 

= — i — T777 rTFn(ni,n 2 )- — — P ril ,n 2 {x,y). (101) 

^—^ n\\n 2 \(N — ni — n 2 )\ ox\ . . . oxNoyi . . . oyN 

Following the same steps, one obtains p 22l p^ and p^. Similarly, one finds that the off diagonal element p 2 ^ is 

1 d 2N 

P23 = T. ~ ~ P23(x,y) 

N\ dx 1 ... dx n dyi . . . dy N 

1 d 2N 

= y2 — j — tttt TrrF 23 (ni,n 2 )- — - — Q ni ,n 2 (x,y) 

n7^ 2 n i- n 2-{ N - ni - n 2 - 1)! oxi .. . dx n dyi . . . dy N 

(102) 

with the polynomial Q ni ,n 2 (x,y) given by 

Q ni .n 2 (x,y) ^^^(wiuD"^^)"^^!^)"^^-^)^-" 1 -" 2 - 1 . (103) 

Using the above equations one can express the density matrix p A in terms of the 2 x N matrix it^fc = ((f>i\(j>B k ) and 
the constants g\ and g 2 . Let us study the particular case where \4>b,) and \4>b 2 ) lie on the plane spanned by \4>\) and 
\<t> 2 )- Then they are related by an SU(2) transformation, i.e. 



> fll ^ _ I cos(0) sin(0)e"A fa 
-&m{0)e- ir i cos(6») J U 2 



(104) 



Making use of (| 102[) one obtains 



p 2 s(t) = — ^= sin(4#)e"'(sin(5i£) sin(V2g 2 t) cos(g 2 t) - s'm(g 2 t) sin(V2gii) cos(g 1 tj) (105) 
2v 2 

which vanishes when cither sin(4#) — or g\ = g 2 . Therefore, in this case, the entanglement transfer scheme works if 
the coupling constants are different. The diagonal elements p\\ and ^44 read 



pn(t) = isin 2 (26»)(cos 2 (%/2 5 it) +cos 2 (V2 52 i)) +cos 2 (20)cos 2 (git)cos 2 ( 52 i) (106) 
,044(4) = cos 2 (26»)sin 2 ( ff it)sin 2 ( 5 2i)- (107) 



From the above equations, one finds that the maximum value of entanglement in AiA 2 is achieved when | sin(4</>)| = 
1. It is interesting to compare this situation with the case when both particles occupy the same state (see (|72|) and 
([75])): if gi 7^ g 2 , one can increase the entanglement transferred to the qubits preparing B in a state of the form (|104[) 
(see FigE]). 



VII. CONCLUSIONS 



We have studied the entanglement induced in a two qubit system as a result of its interaction with a bosonic system. 
The operators coupling each of the qubits to the bosonic system were assumed to commute. As discussed throughout 
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FIG. 5: Negativity as a function of time. The dashed line corresponds to a state of the form (|104|) with 9 = ^ and |J = -j= 
while the solid line corresponds to a 2-particle excitation where each of the particles occupies the state \4>b) = -js{\<t>i) + |</>2)). 



this paper, these interactions appear in the situation when one couples each qubit to a different mode. More precisely, 
we considered operators of the form 

Fi =p l {a i ^a l )a l n + q l {a\a l )cJ j m , i = (1,2) with [ ai ,c^]=0. (108) 

In this case, the mechanism entangling the qubits is analogous to the mechanism responsible for the entanglement 
transfer from two qubit systems to two qubit systems (see Fig l(a)[ ). From section (|V)). we know that a 1-particle 



state being a superposition of the modes \(j>i), \4>2) takes the form of the entangled state 

|*b) = a f (<M |0) = ui |1, 0) + u 2 |0, 1) (109) 

when written in occupation number representation. However, the form of state depends on the interaction 

between the qubits and system B. If one of the qubits interacts with mode \4>b) while the other qubit interacts with 
mode \4>'b) (orthogonal to \4>b)), then the state = a' (</>b) |0) may be written as 

|$ B ) = |1)|0). (110) 

Now, |3>b) has the form of a separable state. It is for this reason that we avoided talking about the entanglement 
between the modes. Instead, we computed the entanglement induced in the two qubit system as a result of the 
interaction with multiparticle systems. For all the N-particles states considered, we found an interaction inducing 
entanglement in the two qubit system. This situation changes dramatically if the bosonic system is in the coherent 
state \z) = e za (4>B)' < -z a(<j> B ) ^ e za {4> B } |q\_ T n f ac ^ thi s state behaves like a separable state for operators of the 
form (|108p . In section (|IV[) . we studied the series expansion of the negativity N{p A ). We computed the first 
nonvanishing contribution to M{p A ) in the case where the operators acting on B were different from those in (|108p . 
We found that when system B is in the particle vacuum state |0), the qubits may become entangled if the interaction 
Hamiltonian contains operators of the form Fj = gi{a{4>i) + (3ia^(tpi)). This type of interactions could be used to 
extract entanglement from a coherent state \z) (entanglement extraction from coherent states has been discussed in 
Q). We leave the this problem for future work. 
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